Large Adiabatic Scalar Perturbations in a Regular Bouncing Universe 
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It has been shown that a contracting universe with a dust-like [w ~ 0) fiuid may provide an almost 
scale invariant spectrum for the gravitational scalar perturbations. As the universe contracts, the 
amplitude of such perturbations are amplified. The gauge invariant variable $ develops a growing 
mode which becomes much larger than the constant one around the bounce phase. The constant 
mode has its amplitude fixed by Cosmic Background Explorer (COBE) normalization, thus the 
amplitude of the growing mode can become much larger than 1. In this paper, we first show that 
this is a general feature of bouncing models, since we expect that general relativity should be vafid 
in all scales away from the bounce. However, in the Newtonian gauge, the variable $ gives the value 
of the metric perturbation 4>, raising doubts on the validity of the linear perturbative regime at the 
bounce. In order to address this issue, we obtain a set of necessary conditions for the perturbative 
series to be valid along the whole history of the model, and we show that there is a gauge in 
which all these conditions are satisfied, for a set of models, if the constant mode is fixed by COBE 
normalization. As a by-product of this analysis, we point out that there are sets of solutions for 
the perturbation variables where some gauge-fixing conditions are not well defined, turning these 
gauges prohibited for those solutions. 

PACS numbers: 98.80.Es, 98.80.-k, 98.80.Jk 



I. INTRODUCTION 

Cosmological models with a contracting phase preced- 
ing a bounce to our present expanding phase have been 
studied as extensions of the standard cosmological model. 
They were analyzed in several contexts, including regu- 
lar and singular bounces [ll-[T^ . In 0, [11, Q , it was 
shown that if the contracting phase is dominated by a 
dust-like fluid and the perturbations are seeded by quan- 
tum vacuum fluctuations, the curvature perturbation C 
develops an almost scale invariant spectrum. 

There are various ways to obtain a bouncing cosmol- 
ogy. For each model there is a specific way to evolve 
the perturbations through the bounce phase. However, 
one can estimate the perturbations in the post-bounce 
phase by imposing some general continuity conditions 
on both background and perturbation variables. It was 
shown in [6] that the curvature perturbation is ampli- 
fied in the contracting phase, and remains constant and 
scale invariant in the expanding phase. Nevertheless, the 
Bardeen [l^ potential $ develops a large growing mode 
in the contracting phase, which is converted entirely into 
a decaying mode in the expanding phase, differently from 
inflation, where the decaying mode is usually small. This 
behavior of the perturbations was also obtained using 
models for the bounce which allows complete calculations 
of the perturbations through the bounce [13 ]. 

Even though the growing mode of 4> couples only with 
the decaying mode in the expanding phase (see for ex- 
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ample j^), its large value at the bounce raises a prob- 
lem concerning the violation of linear perturbation the- 
ory around this time. In a previous work [13,], it was 
explicitly shown, using a specific model, that indeed the 
Bardeen potential grows larger than 1 at the bounce, but 
this mode can be made harmless with a suitable gauge 
choice. Thus remains the question about whether this 
problem is a general feature of bouncing models and, in 
case this is true, whether such gauge choices are still ef- 
ficient to solve this problem in a broader class of models. 
In order to address this point, we show that the ratio be- 
tween the Bardeen potential at the bounce and its con- 
stant value long after the bounce is indeed very large in 
general. Then we obtain a set of necessary conditions, for 
the metric and matter perturbations, that must be satis- 
fied in a valid linear perturbation theory. The first part of 
this set is defined by imposing that the metric perturba- 
tions remain small when compared to their background 
values. The second one comes through the imposition 
that the perturbed Einstein equations remain small when 
compared with the background evolution. As these con- 
ditions are applied to the perturbations, it is necessary 
and sufficient that they should exists at least in one spe- 
cific gauge. This means that it should have at least one 
way of embedding the reference metric in the spacetime 
in which the difference between the two metrics satisfies 
all conditions. Finally, we apply these conditions for the 
perturbations using a specific gauge, and we show that 
they are fully satisfied whenever the constant modes of 
(f> or C in the expanding phase are small. 

This paper is organized as follows. In Sec.[ll|we make a 
brief review of linear perturbation theory around a back- 
ground with homogeneous and isotropic spatial sections. 
One can obtain an almost scale invariant spectrum when 
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the perturbation freezes during a dustlike fluid domina- 
tion. Therefore, in this work, we consider the models in 
which the spectrum of adiabatic perturbations is formed 
within this mechanism. In Sec. lIIIi we discuss the growth 
of adiabatic perturbations in the contracting phase of 
bouncing models with one or more fluids, and the rela- 
tions between the gauge invariant curvature (^) perturba- 
tion and the Bardeen ($) potential. For the bounce cross- 
ing, we review in Sec. IIII Al a nonsingular bounce gener- 
ated by quantum gravity effects in which these two gauge 
invariant variables remain finite and calculable. Then we 
show explicitly that the Bardeen potential becomes very 
large at the bounce. Next, in Sec. IIIIBl we discuss this 
issue for a class of models, and then compare with the re- 
sults obtained in the quantum bounce scenario, arriving 
at the same conclusions. In Sec. IIV Al we obtain the set 
of necessary conditions that scalar adiabatic perturba- 
tions should satisfy in order to keep linear perturbation 
theory valid all along the contracting phase, during the 
bounce, and after it in the expanding phase before recom- 
bination. Then we show in Sec. lIVBi using the uniform 
curvature gauge that all the above-mentioned conditions 
are satisfied along all these phases, except for the bounce 
itself, where in Sec. IIV CI it is shown that one must use 
the synchronous gauge in order to keep all these neces- 
sary conditions satisfied. We end up with discussions and 
conclusions in Sec. IVl 



Defining = p/pc, where pc is the critical density 
today and using the energy conservation equation 
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where Hh = "Hq^ = c/{aoHo) is the comoving Hubble 
radius and Hq is the present value of the Hubble function 
H = a^^da/dt. Note also that H = E/{xKh), where 
E — H/Hq. For a single fluid with w = p/p constant, 
one gets O = ^^qx^'^'^^'^^ and 



/3(I + w) _i 

W, Z = \j — X . 



(7) 



From Eq. ^ we note that Eq. ^ reduces to the equa- 
tion for the perturbations found in [13, for a single 
fluid dominated quantum bounce 



II. LINEAR COSMOLOGICAL 
PERTURBATIONS 

From the Einstcin-Hilbcrt action, one can obtain the 
second order Lagrangian for the Mukhanov-Sasaki (l6j 
variable 
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yielding the equations of motion for their modes with 
wave number fc, 



where 
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In this work, we define n — SttG/c^, T-L = a'/a, a is the 
dimensionless scale factor of the background flat Fried- 
mann model, a prime denotes derivative with respect 
to conformal time dr] — cdt/a, t being cosmic time, 
x = ao/a is the red-shift function, p and p are the to- 
tal energy density and pressure of the matter content of 
the model, respectively, and a subscript denotes the 
present value of the respective quantity; we define z with 
an additional factor Qq ^ compared to that defined in [l6l | , 
Eq. 10.43b. 



< + [wk^ - — ]vk=0, 
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Therefore, we can use this same Eq. ^ to evolve the per- 
turbations in the classical contrac ting phase and through 
this particular quantum bounce (l4l . Il8l | , which we will 
present in the next section. 

Equation ([l} formally looks like a Lagrangian for a free 
scalar field with a time-dependent mass, which can be 
readily quantized. Whenever the potential z" / z becomes 
negligible with respect to the wave number k, vacuum 
initial conditions can be imposed. This happens on the 
onset of inflation or in the far past in the contracting 
phase of bouncing models. 

The Mukhanov-Sasaki variable has the following re- 
lation with the gauge invariant Bardeen potential and 
curvature perturbation [l^, respectively. 



1)2$ = -^/A^lpX^Hz^ (-)' , C = -V^t 
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where Ip = ^jGhlt? is the Planck length and the 
spatial Laplacian. Using Eq. ^ one can also show that 



c = 



2$ 



(10) 



Since the Bardeen potential is a dimensionless quantity, v 
must have dimensions of inverse length. From now on we 
will deal with dimensionless quantities by conveniently 
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multiplying all physical quantities by appropriate powers 
of Rh (e.g. V — > vRh, k — fcR//, ij — >■ rj/Rn, etc.). 
From the operator decomposition,^ we obtain 



Rh V z / 



(11) 



in these situations, which reads 

+ A2{k) ( l-fc2 /'d^e 
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Using the following definitions of power spectrum for the 
Bardeen potential and gauge invariant curvature pertur- 
bation, A| = k^\^k\V{2^^) and = fc3|Cfe|V(2^'), 
respectively, we get 



A<T, = 
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nR^ fc2 
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Ar = ^|-^k'/' 
n Rh 



Note that the amplitude of the perturbations is multi- 
plied by a very small number Ip/Rn — 9.9 x 10~®^ft., 
where h = _ffo/(100 Kms~"'^Mpc~^) is the dimensionless 
Hubble constant. 

The general solution of the mode Eq. (21) can be ex- 
panded in powers of k^ according to the formal solu- 
tion [16,] 
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where we have shown the terms up to order O (/c^) , ry* is 
a convenient choice of the integration constant, related 
to the conformal time where the initial conditions are 
set (in the case of bouncing models we make the choice 
T]* — >■ —oo). The bounce takes place at ry = and any 
function with an over bar refers to its value calculated at 
ry, e.g., f — f{f])- In Ec^. (fT3)) . the coefRcients Ai and A2 
are two constants depending only on the wave-number k 
through the initial conditions. Once the solutions freeze 
(Cj/c^ ^ z"/z), i.e., when the mode is below its poten- 
tial, the superhorizon solutions above can be used. Since 
we are interested in what happens with the amplitude 
of the spectrum after it is formed in a dust-dominated 
evolution, it is enough to analyze the superhorizon so- 
lutions. For long wavelengths of cosmological relevance, 
this happens during the bounce, of course, and around 
our expanding epoch. Under these conditions, we can 
use Eq. (|T3l) up to second order in k^ to calculate {v/ z)' 



^ The creation/annihilation operator decomposition of <& involves 
terms like ^i^a^^d^k. Since the operator <I> is dimensionless and 
aj^ has unity of length L^^'^ , <J>fe also have dimension of length 
L3/2. 



III. THE PROBLEM 

In this section, we will calculate A$ for long wave- 
lengths using Eq. (HH) after the decaying mode be- 
comes negligible in the expanding phase A^ and near 
the bounce A|,. We will show exactly for the quantum 
bounce of Refs. and using general arguments for 

general bounces that A^ is many orders of magnitude 
larger than A^. Nevertheless, in the next section, we 
will show how the linear perturbation theory is still reli- 
able near the bounce in spite of the problem, which will 
be described in the sequel. 

Let us evaluate A|, and A^ . In the case of A|, , we can 
see from Eq. (fTl)) that the term multiplying Ai (fc) is the 
decaying mode of the contracting phase, which goes to 
zero at the bounce. Hence, when the solution gets close 
enough to the bounce, the important contribution for A| 
is 
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For A^, the term multiplying Ai{k) is the growing mode 
of the expanding phase. However, the last term which 
multiplies A2{k) can be written as 
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where we defined the constant B = dryz"^. We also 
have discarded the last term in the sum because it cor- 
responds to a decaying mode in the expanding phase, 
and by the positivity of the integrand, it will always be 
smaller than B. 

The presence of i?, as we will see, makes the term 
multiplying A2{k) much more important than the one 
multiplying Ai{k) in the evaluation of A^. Hence, we 
get 
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when the decaying mode becomes negligible. 

The variable C, has a much simpler evolution, in the 
contracting phase, the term 
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grows as the perturbations approach the bounce. In the 
expanding phase the integral can be spht as in Eq. (ITSl) . 
where the first term is constant and the second becomes 
the decaying mode. Unhke this decaying mode around 
the bounce has the same order of magnitude than its con- 
stant mode. Thus, after the bounce the power spectrum 
of C is 



l^k'/^ \A,ik)B\ 



(18) 



Generally, as we will show in the next sections. 



df]clz' = O (1) . 



Therefore, using Eqs. P7)) and (|18p . we note that in the 
expanding phase sa A^. This shows that the power 
spectrum of C is enough to assess the spectrum of $ in 
what concerns its constant mode. However, it is insen- 
sitive to the large growing/decaying mode developed by 
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We will now evaluate A^ and A|, for a particular quan- 
tum bounce [13, [3 and then for a more general case. 



A. The example of the quantum bounce 

References [3, [l^ [2^ show the canonical quantiza- 
tion of a minisuperspace cosmological model describing 
a perfect fluid with p = Wqp on a Friedmann geome- 
try with flat spacelike hyper-surfaces. The corresponding 
Wheeler-DeWitt equation reads 
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I 9^vl/(0)(a,r) 
4 dx^ 



(19) 



where x = a^^^""'''-'/^2/(3(l-Wq)). This is just the time- 



= a3(l-«.,)/2r 

reversed Schrodinger equation for a one-dimensional free 
particle constrained to the positive axis. 



After imposing a Gaussian initial wave function, the wave solution of this equation for all times in terms of a reads 
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where T is related to conformal time through d?7 — [a{T) 
width of the initial Gaussian. 



^ AT and T;, is an arbitrary constant related to the 



Because of the chosen factor ordering, the probabil- 
ity density piii(a,r) has a non trivial measure and it is 
given by p*(a,T) = a(i-3«.<,)/2 l^^^^^^ 2.)|2_ continu- 
ity equation coming from Eq. (jll 



\ reads 



dT 



d_ 

da 



2 'd^ 



0, 



where S denotes the imaginary phase of the wave func- 
tion ^fnv This implies in the de Broglie-Bohm interpre- 
tation mHH that 



da 
df 



2 9^' 



(22) 



in accordance with the classical relations da/dT — 
{a,H} = -a(3'^''-i)Pa/2 and Pa = dS/da. 

Inserting the phase of Eq. (pOt into Eq. ([^ , we obtain 
the Bohmian quantum trajectory for the scale factor: 



a{T) = flfc 



1+ ^ 



TV 



(23) 



Note that this solution has no singularities and tends to 
the classical solution when T — >■ ±c». Solution (l23l) ca n 
be obtained for other initial wave functions (see Ref. jl9|). 



Changing to cosmic time dt = a'^^'^dT, we obtain the 
Hubble function 

1 da _ 2ra-3""' 
adt^ " 3(T2 + T2)(1-u,,)- 



(21) Solving Eq. ^ for T, 



r = ±TM — 




3(1— uj,) 



- 1, 



we obtain 

4„-6'"<i„-3(i-'",) / 6 



9{l-w,fT^ 



(24) 



X 



where Xh = ao/af,. 

Equation (j24p is equivalent to the Friedmann equation 



(25) 



with the additional term —HQ^qox . The quantity = 
flqox^ is an effective density related to the quantum evo- 
lution. Thus, the evolution of the scale factor is equiva- 
lent to the one obtained by adding to the matter content. 
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in the classical Friedmann equation, a stiff negative en- 
ergy fluid, i.e., p + ps, where ps = -VLqopcX^'^'^+'"''\ 
Note, however, that this is a quantum effect and there- 
fore there is no perturbation associated to this effec- 
tive fluid. Note also that away from the bounce phase 
{H^VLqax^ < Kc2p/3) we obtain H"^ = i/|f7^o2;^^^+'"'\ 
where 



9„-3'^9^-3(l-™g)/2 
3(1 - Wq)Tb 



(26) 



Comparing Eqs. ([M)) and (pS)) . we obtain the impor- 
tant relation for the dimensionless T;,, 



Th = 



2ar 



3(1 - Wq)y/Tr^ 



(27) 



In this bouncing model there is one single fluid. Hence, 
Eq. (0 holds. Using Eqs. and ([271), we obtain 
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and 
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where flyj — flowX^'"^'^^''''; in the last equation we have 
neglected the evaluation of the primitive at ry = 0. Note 
that in this simple model one must have Wq ~ in order 
to obtain an almost scale invariant spectrum of pertur- 
bations. Hence, as Xb — ao/at » 10^° (we expect that 
the bounce occurs much earlier then nucleosynthesis), the 
quantity B is indeed a large number. 
From Eqs. (fT5|) and (TTfl) we obtain 
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and 
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AlpV2^Wqk^/^ \A2{k)\ 
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As expected, does not depend on time. The ratio 
between these two quantities is 

^OiL.3(l - Wq){5 + 3Wq) l+3w. 



A'' 
A° 



Ml 



(32) 



As the Cosmic Microwave Background CMB observa- 
tions [Hi require ^ 10"^ and xt > 10^°, then 
A|, :S> 1, which turns questionable the validity of lin- 
ear perturbation theory at the bounce. We will now see 
that this issue is also present in a much larger class of 
bouncing models. 



B. More general bounces 

The general solution for the Mukhanov-Sasaki vari- 
able, Eq. (TT^ . for the adiabatic perturbations, is valid in 
the contracting and expanding phases when the dynam- 
ics are given by general relativity (GR) and also through 
the bounce itself, in the case of the quantum bounce dis- 
cussed above. For general bounces, one is not sure that 
the solution in Eq. ([T^ is valid through the bounce due 
to not having any particular analytic solution in order 
to evaluate it away from the bounce, as we did in the 
last subsection. However, if the bounce is short enough 
an estimate of Eq. (1131) away from the bounce, when GR 
is valid, will be sufficient to evaluate the orders of mag- 
nitude of Eqs. (ITSl) and ([T7| as long as a short bounce 
does not change these figures too much due to the ex- 
pected continuity of perturbations through it (see [l^ 
for a general discussion about matching conditions us- 
ing the continuity of the perturbations and ^] for its use 
in this context). However, one could have a bouncing 
model with a long-time scale. Thus, in our analysis we 
are assuming that the characteristic time of the bouncing 
model is small enough (usually of the order of Ip) that 
we can ignore this phase. 

In order to estimate (flSl) and p7)) we must evaluate 
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For i?, we first divide the integral in the pre- and post- 
bounce branches 



°%.r'% (33) 

-oo ^ Jo 

In each branch the scale factor and, consequently, the 
red-shift variable x can be used as a time variable. Per- 
forming this transformation, we have 

r° dr? _ n dx n dx 



dr] 

T2 



E-z'^ 
dx 



\E-\z-- 

dx 

E+z^' 



(34) 
(35) 



where E~ is the dimensionless Hubble function H/Hq 
during the contracting phase and, therefore, a negative 
quantity; E~^ is the dimensionless Hubble function during 
the expanding phase. We have assumed that a;(±oo) = 
but this is not necessarily true. However, one can always 
assume that x{±oo) / Xb ^ 1, which is sufficient to esti- 
mate the integral. Also, an asymmetric bounce will not, 
in general, change too much the orders of magnitude we 
will evaluate. Hence, we will assume for simplicity that 
the bounce is symmetric 



dx 



E+z^ 

Now we divide the above integral in two domains, 

da; 



r" dx 


r" dx 




L Ez^~ 


Jo Ez^^ 


I 

J X 



(37) 



6 



where Xc is the value of the red-shift function where the 
new physics of the bounce begin to be relevant {xc could 
be defined as the value of x in which d^a/Afi — 0, the 
transition from the decelerating behavior typical of GR 
to the accelerating phase of the bounce. Note that the 
second portion of the sum above cannot be written if the 
solution ([T3| is not valid through the bounce. However, 
as Xc is generally of the same order of magnitude as Xf,, 
then Xc S> 1 (see the quantum bounce example, where 
Xc = [(1 + 3w5)/4]i/[3(i-«'i)]a;^ ~ xb for < Wq < 1). 
In this case, the interval Xc < x < Xb '^s irrelevant when 
compared to < a; < Xc and, therefore. 



dx 



dx 



(38) 



In the GR domain, we can use Eqs. Q and ([5]) to 
obtain 



2cix 



.2^3 



iEil+p/p)' 



(39) 



As X increases, the quantities above become dominated 
by the fluid with largest value of Pq/ Pq — Wq, and in this 
phase the integrand is dominated by the term a:'^'^^"™'')/^ 
since in one fluid domination and p/ p become constant 
and E cx x^^^+^'i)/'^. For simplicity, we are assuming 
that the fluid which dominates in this epoch has constant 
equation of state. As x varies several orders of magnitude 
in the integration interval, the value of the integral is 
dominated by the integrand near instant Xc, where we 
are assuming that the fluid with equation of state Wq 
dominates. Using this feature, we show in the Appendix 
that this integral is approximated by Eq. (IA.4I) . i.e.. 



dx 

'e^ 



3(1 - Wq) E{x,;)z{xcy 



(40) 



where, as justified before, we are taking Oc « flf,, and we 
have used Eq. dSH]) for the last equality. The integrand 
in the last integral is an increasing function of a, hence, 
using the mean value theorem, we get 



da 



-? -2 

a?E 



aoc^z 



a^E 



(a - Oc) 



< 



E 



(43) 



where in the last approximation we used a 3> flc and 

Now inserting Eqs. (gT]) and dU) into Eq. ([T7]), assum- 
ing for simplicity that near the bounce there is domina- 
tion of one fluid, we get 

^ [2 AlpWqky^ \A2ik)\ 3(l-^,)/2 , . 

V TT 3(1 - WqjRHV^wO 

which, apart from numerical factors of order unity, coin- 
cides with Eq. dnH). 

For the general one has 



A'' 



and the ratio is 



A'' 
AO 



> 



^Ow3{l — Wq) l+3u,, 
4(^fc)2 



(45) 



(46) 



Again, this ratio is proportional to x[^^^™''\ and for 
any fluid with Wq > —1/3 (assuming COBE normaliza- 
tion) one gets a very large amplitude during the bounce. 
As we have shown, this large ratio A^/A^ is a general 
feature of a long contracting phase and, therefore, it will 
be present in any reasonable bouncing model where the 
matter content satisfies —l/3<p/p<l. 



As Xc ~ Xb, we obtain 



B : 



2xb 



3(1 - Wq) E{xb)z^{xb) 



where it must be understood that, although evaluated at 
Xb, the functions E{x) — \/fl{x) and z^{x) in Eq. (|1T|) 
are the usual GR expressions for them, which are valid 
just before the bounce. 

Note that if the integral J^"" dx/{Ez^) makes sense dur- 
ing the bounce, it is a positive quantity that is being 
neglected in the evaluation of B and, hence, Eq. 
continues to be valid. Note also that although E — at 
the bounce, this integral converges for regular bounces 
(see the quantum bounce above). 

For the second integral one has 



1 — 2-2 

drjCgZ 



da 



-1 -2 

a^CgZ 
a?E 



da 



-2 -2 

a?E 



da 



3(1+P/P) 
2aQE 



IV. THE SOLUTION 

(41) 

As we have seen, the gauge invariant Bardeen poten- 
tial $ may grow in the bounce because what would be 
the decaying mode in the expanding phase is the grow- 
ing mode in the contracting phase. This mode can be 
very large around the bounce if the contraction is huge. 
Then one could put into question the validity of linear 
perturbation theory at the bounce, which compromises 
all calculations of cosmological perturbations in bounc- 
ing models. However, the definition of a gauge invariant 
quantity is not unique. One simple reason for this is 
that one can multiply any gauge invariant quantity by a 
background function and it continues to be a gauge in- 
variant quantity. For instance, if one defines the gauge 
invariant function (a/ao)"$ with a positive, which co- 
incides with the Bardeen potential today, it is trivial to 
find a power a, where this gauge invariant function is 
small at the bounce. Hence, what one has to do is to 
(42) look at Einstein's equations for the perturbations them- 
selves and see if the linear theory makes sense, at least in 
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some gauge, during the whole history of the model before 
the usual epoch where nonlinearities become important. 
Note that it is not necessary that the theory makes sense 
in all gauges; a valid gauge transformation which relate 
different gauge choices at some phase in the cosmological 
evolution may not exist. 



A. Conditions for linearity 

Let us concentrate on the scalar perturbations. The 
geometry of spacetime is given by 



.9p 



- a'-°^ + h 



(47) 



where g^J represents the homogeneous and isotropic cos- 
mological background 

ds2 = g(°Jdx''dx'' = -c^dt^ + a^-f^jdx'dx^ , (48) 

where 7y is the metric of the maximally symmetric 
spatial hyper-surfaces with normalized scalar curvature 
K = 0, ±1, and /i^^ represents linear scalar perturba- 
tions around it, which we decompose into 



hoo = 20, 
hoi = -aDiB, 



(49) 



where Di is the covariant derivative with respect to "fij. 
Hence, the first conditions are 



< 1, 6< 1, 1, V< 1- 



(50) 



Computing the perturbed Einstein equations in the 
background Gaussian coordinate system, one gets 



2ese sn , 

SGo = — = -Kdp kSTq , 

= Id, ( Se - ^5<j - ^5<j 



(51) 
(52) 



A[k(p + p)V] = K<5I^, 







SGi 



-2( — + 659 + 0— -Aa^ + ?h (53) 
a a 4 



— 3kSp — KSTi' 



[i:~(t> —5(T 



i2 



(54) 



-(5n = n5Tj\ 



where in the fourth equation i ^ j, = DiD^, Sp, 
and Sp are the perturbed energy density and pressure, 
respectively; 511 is the anisotropic pressure, which we 
will consider to be null, and V is the perturbed velocity 
field potential. 

In these equations, the quantities 5Q, Sa, STZ and a* 
are the perturbed expansion rate, shear, curvature scalar. 



and worldline acceleration with respect to the constant 
cosmic time hyper-surfaces. The background expansion 
rate is simply & — n'^-fi — 3H/c — 3H/a, where is 
the normal of the maximally symmetric spacelike hyper- 
surfaces and ■ represent the covariant derivative compat- 
ible with the background metric. 

These quantities are related to the metric perturba- 
tions by 



Sa = 

sn = 



-D^{£' -B) 
4 



(55) 
(56) 

(57) 
(58) 



The gauge invariant variables are defined by the following 
combinations of the perturbation variables: 



= + 



a 



(59) 



We have to verify whether the perturbed Einstein 
equations remain small when compared with the back- 
ground Einstein equations, where the non-null back- 
ground Einstein's tensor components are 



n _ 

Lrn — — 



G,' = -7f 



82 
3 



3K 



K) 



K 



26' 
3a 



a 
3 



2%' + 



(60) 



(61) 



For Eq. dCT]), the term 959 is related to 9^ in G^P and 
5TI is related to K/a^ in Gq^. Hence, one must have 



1591 « 191 



a59 



« 1, 



and 



{D^ + ■iK)ip 



K 



< 1. 



(62) 



(63) 



In the case where K 
Go° and Gi", yielding 



«1, 



0, we have to compare 5TZ with 



2W + 



< 1. 



(64) 



Now we have to establish the conditions on 5a and a*, 
which are null in the background. We will use Eqs. ([55)1 
and which come from perturbing Gj*. These com- 
ponents of the Einstein tensor contain, in the back- 
ground, 9'/a and 92. The first one originates 5a' /a 
and Dia^ in these equations, while the second originates 
the term Q5a in Eq. (|55)1 . Thus, we obtain the following 
conditions: 



DHa 



an 



« 1, 



Z?2 



n'-w 



< 1. 



(65) 
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There are no further independent conditions on the geo- 
metric perturbations. Note also that each condition was 
obtained comparing the background and perturbed val- 
ues of each term of the perturbed Einstein equations. 
This means that these conditions are stronger than just 
comparing the perturbed equations with the complete 
Einstein tensor. 
The conditions 



« 1, 



< 1, 



when in terms of geometric perturbations using Eqs. (j5ip 
and ((53)) . reduce to the conditions already obtained for 
those quantities. The velocity field potential V appears 
only in the JGi" projections, which are null in the back- 
ground. To deal with this variable we can calculate 
the second-order correction on the trace of the energy- 
momentum tensor, which contains a term proportional to 
DiV D^V . Such a term also appears when one calcu- 
lates the normalization of u^jW^ = —l+a~^DiVD'^V up to 
second-order terms, where i;^ represents the eigenvector 
of the energy-momentum tensor; we have already con- 
sidered the other quadratic corrections in this expression 
much smaller than 1. Therefore, applying the condition 
a"^|AV| < 1 on Eq. ([52]), we obtain 



aD.Se 



W + K-W 



« 1, 



a{W + K-n') 



< 1. (66) 



These conditions are weaker than Eqs. and ([55)) 

around the bounce phase, first due to the additional fac- 
tor 'H~^ (X rj, which goes to zero near the bounce, and 
also because the spatial derivative brings down higher- 
order correction terms when acting in the superhorizon 
solutions (Eq. (IHl)). 

The Weyl tensor C^ua^ in a Friedmann metric is 
null. Therefore, its perturbation will be gauge invari- 
ant (see [13|)- One could expect to get a gauge invari- 
ant condition for the perturbative series by comparing 
the perturbation of the Weyl tensor with the background 
Riemann tensor. The non-null components of the per- 
turbed Weyl tensor are given by its electric part 

<SC,o^° = a-^ (^D,D^ - ^D^^ ($ + vj/), 

while the background Ricci tensor components are 



2K 



6' 62 



, 2K + 'H' + 2n^ 



Comparing these components yields the constraint 



2K + n' + 2-H2 



«1, 



which is satisfied whenever conditions given in Eqs. (|63f - 
[65)) hold. Hence, the Weyl tensor provides a gauge in- 
variant constraint, which is necessary but not sufficient 



since it alone does not imply Eqs. (|63l - [65)) . Additionally, 
this condition, when applied to the perturbations near 
the bounce phase is much weaker than that of Eq. ()50l) , 
which is necessary to define the perturbations of the in- 
verse metric. 



B. The gauge choice solution 

The evolution of the gauge invariant Bardeen potential 
in the classical GR phase around the bounce is described 
in Sec. IIII B[ where it shown that $ grows larger than 1 
at this phase. In the Newtonian gauge {B = = £ and 
consequently 6a = 0), $ = 0, therefore condition (|5(I)) is 
not satisfied and linear perturbation theory breaks down 
in this gauge. 

However, choosing a gauge with constant curvature 
{6TZ — 0), we avoid the problem described above. In this 
gauge one sets ip — B — and, in order to completely 
fix the gauge, we impose £{rii) = 0, where r]i is some 
particular convenient choice of conformal time. One has 



Sa = 



IT' 



2 2 2/- 



d^$ 



(67) 



where we used $ = \I> {5Ii = 0) and Eq. (dU)) to obtain 
the expressions above. 

In this gauge the perturbation (/) has a different behav- 
ior. At any instant in which a single fluid dominates one 
has 



'S{l + w) 



c, 



where we used Eq. ([7]). Therefore, in this gauge cf) fol- 
lows the evolution of ( instead of $. As we discussed in 
Sec. lIIIi C, grows in the contracting phase until it attains, 
near the bounce, an amplitude approximately equal to 
the constant mode of In the expanding phase, C also 
has a decaying mode, but in this case this mode is always 
smaller than the constant one. Hence, <C 1 is satisfied 
in this gauge. 

Starting the calculations in the constant curvature 
gauge, one can see that near the bounce scale the gauge- 
fixing condition for the Newtonian gauge is not well- 
defined. In the new gauge the value of 4> would change 
as 3{1 + w)/2(^ $. However, as we showed near the 
bounce $ 3> C- This would imply a non valid transfor- 
mation 



<^N = 0CCG + 'I' 



3(1 



-c 



where 0n represents the metric perturbation in the New- 
tonian gauge and 0ccG in the constant curvature gauge. 
The relation between the perturbation £ and $ is 



\£\ = 



{xi - x) 



< 



E^Xb) 



9 



where we used the mean value theorem. Note from 
Eq. (gS]) that near the bounce ^{xb) oc a;(5+3"'g)/2 a^d 
thus x'^/E'^ cx in this gauge £ oc x^^^"'""'/^ 

has the same growth factor as (f>. With the results above 
and noting that in this gauge tp — 13 — 0, the conditions 
given in Eq. (jSOp are verified. The other conditions given 
in Eqs. ((62H65|) can be trivially verified. 



gauge fixing condition. However this approach alone is 
not enough to evaluate the validity of the linear approx- 
imation. Therefore, one must always check if there is a 
gauge in which the perturbation series is valid. 
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C. Gauge choice and the bounce phase 

The discussion above shows that for a variety of models 
the evolution of the perturbations near a bounce phase 
is well-behaved in what concerns the perturbative series. 
However, at the exact moment when the bounce occurs, 
other problems can arise. The gauge invariant variable C 
is related to the metric perturbations through (Eq. EZl) , 



a^K(p -I- p) 
2^2 



Hence, at the bounce the Hubble function goes to zero 
and, therefore, the perturbation cji diverges since C stays 
constant at the bounce (see Eq. [T3|) . However, in the 
synchronous gauge one has 



n 



dr] 



2^2 



(68) 



and one can show that, using the variable /i defined as 
, one has H oc ^ near the bounce (this is 



XbC 



the case when the bounce is caused by a negative fac- 
tor in i?2). Thus, the integral above is proportional to 
/J, / d^/i^2^ and the perturbations are well-behaved at the 
bounce. Hence, the perturbations are always finite and 
small, as the constant mode of $. 



V. CONCLUSIONS 

We have shown in this paper that, for adiabatic pertur- 
bations, the Bardeen potential in the contracting phase 
of bouncing models can generally become very large, but 
this fact does not invalidate linear perturbation theory 
around the bounce. We established necessary conditions 
for the validity of linear perturbation theory on Fried- 
mann backgrounds, and we have shown that there is a 
gauge choice, for a large class of bouncing models, where 
these conditions are satisfied. In fact, there are some 
gauges that are ill-defined close to the bounce because 
the gauge transformations relating them to some well- 
behaved gauge are not valid. 

In conclusion, the program of describing the evolu- 
tion of linear primordial perturbation in bouncing mod- 
els is well-defined. However, one must take care with the 
gauge which will be chosen while performing calculations, 
since some gauge fixing conditions are not well-behaved 
in these scenarios. In general, the gauge invariant ap- 
proach is more appropriated. It does not depend on any 
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Appendix: Asymptotic Series Approximation 

The integral ([55]) can be rewritten as 

(A.l) 



dx 



— - / dz.5(i.)e-. 



where we have defined v — ln(2;), s = 3(1 — Wq)/2, = 
ln(a;c), and g[y) = Xj {x'^^-'^'^'^'^l'^Ez'^). This divides the 
integrand in the dominant term e"'" ~ ^3(i-u)<j)/2 
the controlled function giy). Using Eq. (HH]), it can be 
written as 



1 



2c2^3(l+«,,)/2 



2,(l-3«;,)/2^^2 ZEiX+pIp) 



(A.2) 



Hence, 



lim g{v) = O (I) , lim g{v) = 0. 



The first limit comes from the fact that E oc x-^^^+^'^Z^ 
for X — >■ Xc. This could also be seen in the situation where 
the matter content, besides the fluid with Pq/ Pq = Wq, 
is given by a collection of n other constant Wi fluids, 
yielding 



3{1+P/P) 



1=1 



-1/2 



As Wq > Wi, when in the domain in which x 3> f , gii") is 
dominated by a constant value of order one. 

The second limit corresponds to x — or a — >■ oo, and 
one expects that the universe was dominated by the fluid 
with wi since we ordered the fluids imposing Wi < Wi+i, 
giving g{iy) cx 2:3(io,-«.i)/2 _^ q 

With these results, we can integrate Eq. (jA.ip by parts 
to obtain 



dx e'^"" , , 
Ez^ = — 



dv 



e^" dg{v) 
s du 



(A.3) 



Note that 



lim x''!^ = 0(1), lim x'--^=0. 
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where r = 3(wg — Wn)- Thus, we can again spht the 
integrand in a controUed function gi = x^dg/dv times 
g(s-r)i/^ and integrate by parts obtaining 



-9{vc) 



e'^"" dg 
s(s — r) dv 



dv 



The magnitude of the second term has an additional fac- 
tor of x~^ when compared with the first. Repeating the 
process of factoring the largest growing term and inte- 
grating by parts, we obtain an asymptotic series approx- 
imation for this integral. For a complete discussion about 



this method of approximating integrals, see [28| . 

As we have shown for the first term, each subsequent 
term will by multiplied by an additional factor of Xc to 
a negative power. Therefore, to estimate the order of 
magnitude of the integral, it is sufficient to keep only the 
first term, 



dx 



Sz2 3(1 - Wq) E{x,)z{xcy 



(A.4) 



It is worth noting that the result above is exact if the 
matter content consists in just a single fluid with constant 
equation of state. 
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